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Abstract. I present the first results of a search for clusters of galaxies in Chandra ACIS 
' pointed observations at high galactic latitude with exposure times larger than 10 ks. The 



survey is being carried out using the Voronoi Tessellation and Percolation technique, which 



\l ' is particularly suited for the detection and accurate quantification of extended and/or low 

o ; 

^ surface brightness emission in X-ray imaging observations. A new catalogue of 36 cluster 



candidates has been created from 5.55 deg^ of surveyed area. Five of these candidates have 
already been associated to visible enhancements of the projected galaxy distribution in low 
deepness DSS-II fields and are probably low-to moderate redshift systems. Three of the 
candidates have been identified in previous ROSAT-based surveys. I show that a significa- 



C/3 

. tive fraction (30-40%) of the candidate clusters are probably intermediate to high redshift 

systems. In this paper I publish the catalogue of these first candidate clusters. I also derive 
•i-H . 

■ the number counts of clusters and compare it with the results of deep ROSAT-based cluster 

■ surveys. 
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1. Introduction 

Clusters of galaxies are very important objects in astrophysics. Observations of such systems at 
different redshifts can be used, for instance, to explore galaxy evolution in very dense environ- 
ments (e.g. Dressier et al. 1997; Jorgensen et al. 1997; Kelson et al. 1997; Ellingson et al. 2001; 
Nelson et al. 2001; van Dokkum & Franx 2001). Clusters are also crucial probes of mass dis- 
tribution at intermediate and large scales (e.g. Nichol et al. 1992; Collins et al. 2000; Schuecker 
et al. 2001; Gonzalez et al 2002). Moreover, in the scenario of hierarchical structure formation, 
the estimate of cluster abundance as a function of their mass and redshift plays a fundamental 
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role in constraining the cosmological parameters and erg (e.g. Baiicall & Cen 1993; Girardi 
et al. 1998; Borgani et al. 2001; Reiprich & Bohringer 2002). Following a different approach, 
observations of the Sunyaev-Zeldovich effect (Sunyaev & Zeldovich 1972) in a large sample of 
distant clusters can also be used for a direct measurement of their distances, and can thus provide 
estimates of the Hubble's constant (Mauskopf et al. 2000; Reese et al. 2000). 
In any case, all these studies require large samples of clusters in a wide range of redshift, from 
z = to z = 1 and beyond. In particular, nowadays much interest is given to the search for 
medium-high redshift clusters, whose physical and statistical properties need a big effort of in- 
vestigation. Up until the present, the largest samples of distant clusters have resulted from op- 
tical/NIR and X-ray sky surveys. Optical/NIR cluster surveys consist in searching for enhance- 
ments in the surface density of galaxies (e.g. Postman et al. 1996; Nonino et al. 1999). This 
classical technique suffers seriously from projection effects, sometimes leading to false identi- 
fications, particularly for poor clusters at high redshift, where the contrast with the background 
galaxy surface density is low. Another problem of these surveys is that the cluster galaxy rich- 
ness has a poor correlation with mass (e.g. Borgani & Guzzo 2001). This does not assure that 
the selected sample contains all clusters within a given volume and above a given mass, since 
the cluster selection function is only loosely related to this fundamental parameter, which is what 
model predictions are based upon. Therefore it is difficult to use optically/NIR selected clusters 
for cosmological studies. The advantage of optical surveys is the availability of a large number 
of ground telescopes with high quantum efficiency detectors. This makes completing wide area 
surveys an easy and relatively quick task. In fact, recent deep optical/NIR surveys (e.g. the Las 
Campanas Distant Cluster Survey) are now reaching areas in excess of 100 deg^ (Gonzalez et al. 
2001; Nelson et al. 2002) enabling the discovery of the richest clusters at z ~ 1 and the study of 
the evolution of the moments of the cluster galaxy distribution. 

Instead, X-ray cluster surveys consist in detecting the X-ray emission from the hot gas which 
fills the space between cluster galaxies. Even if large area X-ray surveys are time intensive, their 
big advantage is that clusters look sharper in the X-ray than in the optical sky. In fact, while the 
optical emission is approximately only linearly dependent on the number of cluster galaxies, the 
X-ray emission is proportional to the square of the local gas density. Moreover, the detection of 
the intracluster medium is an unambigous indicator of a real cluster. For this reason the prob- 
lem of spurious detection is much less severe. But the best feature of X-ray surveys is that their 
selection function can be determined to high accuracy knowing the properties of the telescope 
used. Moreover, since X-ray luminosity is a good indicator of the total cluster mass (Reiprich & 
Bohringer 1999), we can safely make cosmological predictions for clusters of a given mass and 
sensibly compare them with observed clusters of a given luminosity. In synthesis, X-ray surveys 
are much easier to use than optical/NIR ones as cosmological probes to address global cluster 
properties and not only for studying single, yet interesting, objects. 

In the past decade these good qualities encouraged the compilation of many X-ray cluster surveys 
mainly based on ROSAT data (e.g. Rosati et al. 1995; 1998; Scharf et al. 1997; Vikhlinin et al. 
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1998). All these surveys gave important contributions in determining the X-ray cluster number 
density and luminosity function at lower and lower flux limits. 

In this paper, I present the first results of a search for galaxy clusters in deep (exposure time T>10 
ks) Chandra ACIS pointed observations at high Galactic latitude. The Chandra's high spatial res- 
olution is the strong point of a survey based on this telescope, since it is essential to efficiently 
separate between pointlike and extended sources and, in particular, to reduce the problem of the 
contamination of distant clusters by active galactic nuclei. The survey is being carried out using 
the Voronoi Tessellation and Percolation technique (VTP, Ebeling & Wiedenmann 1993). This 
approach is completely non parametric, i.e. it is equally sensitive to spherically simmetric struc- 
tures and to irregular ones, and shows its effectiveness to detect low surface brightness sources 
with greater sensitivity when compared with standard detection methods (Scharf et al. 1997). 
This paper is arranged as follows: in Sect. 2 I present a brief overview of the Chandra selected 
pointings and the data reduction; in Sect. 3 I describe the source detection method and the steps of 
the construction of the catalogue. In Sect. 4 I present the calibration of the survey sky coverage. 
Sect. 5 presents an estimate of the redshift distribution of the detected candidates and a measure 
of the sky volume probed by the survey. In Sect. 6 and 7, 1 present the log A^-log S relation and 
the catalogue of detected candidates. Finally, in Sect. 8, I summarize the main results and give 
my conclusions. 

Computations throughout the paper use Qm = 0.3 and Qa = 0.7. The Hubble constant is Hq = 75 
km s"' Mpc"', unless otherwise stated. All X-ray luminosities and fluxes are reported in the 0.5- 
2.0 keV energy band. 

2. Observations and data analysis 

The data used in this paper are 81 Chandra ACIS pointed observations downloaded from the 
public archive. Selected observations have exposure times larger than 10 ks (see Fig.^ and ab- 
solute galactic latitude larger than 20 deg. This sky-selection criterium minimizes the effects of 
absorption of extragalactic X-rays from the galactic neutral hydrogen concentrated around the 
Milky Way plane. In general, the galactic hydrogen column density along the line of sight of 
the selected pointings is nn < 10^' cm"^ (see Fig.|2ji. I preferentially avoid pointings whose 
targets are bright extended sources (typically nearby clusters and galaxies), since they dominate 
the field of view and make it impossible the detection of other objects. This also ensures that the 
cluster selection function is not significantly affected by the known angular correlation between 
galaxy clusters. In Fig.|3j the positions of the 81 selected Chandra fields are plotted in Galactic 
Coordinates. The size of the points is proportional to the exposure times, which ranges from 1 1 
to 168 ks. 

The Chandra ACIS is a scientific instrument made of two blocks. One, ACIS-Imager (ACIS- 
I), is a square matrix of 4 1024x1024 pixels front-illuminated (FI) CCDs (chips 10, II, 12 and 
13) covering a field of view of about 17 by 17 arcmin and is used for imaging observations. 
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Fig. 1. Histogram of the exposure times for the 81 selected Chandra ACIS pointings. 



The other one, ACIS-Spectrometer (ACIS-S), is an array of 6 1024x1024 pixels CCDs (4 front- 
illuminated, 2 back-illuminated (BI)) covering a field of view of about 50 by 8 arcmin and is 
used both for imaging observations and, mainly, for spectroscopical observations with gratings. 
In general, during an observation ACIS can operate simultaneously both with ACIS-I CCDs and 
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Fig. 2. Histogram of the galactic neutral hydrogen column density along the line of sight of the 
81 selected Chandra ACIS pointings. 
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with ACIS-S CCDs, for a maximum of 6 active CCDs. I consider here both images taken by 
ACIS-I chips (chips 10, II, 12 and 13), and ACIS-S imaging observations (aimpoint on S3) taken 
by chips S2, S3, 12 and 13. In ACIS-S observations I take into consideration only the FI chips 
(S2, 12 and 13). I choose to search for X-ray sources only in these chips since they manifest a 
similar photon background. This is important for a correct use of the source detection algorithm 
(see Sect.E}- Moreover, in the selected ACIS-S pointings the S3 chip is often completely filled 
by target photons and, therefore, unusable for the search of other sources. 
Before Chandra observations be used to search for sources, they need some reduction. I per- 
form the data reduction using the dedicated software CIAO (Chandra Interactive Analysis of 
Observations, Fruscione & Siemiginowska 2000). For each pointing, first I filter the preprocessed 
event files selecting only events with a status equal to zero and with ASCA grades 0,2,3,4 and 
6. This particular grade selection appears to optimize the signal-to-background ratio of sources 
(Chandra X-ray Center Team, 2001). Then, I clean bad offsets intervals and examine the data on 
a chip by chip basis. Since for each photon event the energy is known, I select all events with 
energy between 0.3 and 10 keV. This is the band in which the energy scale of ACIS CCDs is cal- 
ibrated. For each chip I filter out bad columns and pixels and remove times when the count rate 
exceeded three standard deviations from the mean count rate per 3.3 s frame time intervals. In 
fact, the background is generally quite stable during an observation. Nevertheless, occasionally 
there are significant variations, especially in the BI chips but also in FI ones. Since for each pho- 
ton event the occurance time is known, I am able to filter out high background intervals. I then 
clean each chip for flickering pixels (pixels showing unstable states of dark current signature), 
i.e. I remove times where a pixel had events in two sequential 3.3 s intervals. At the end, I remain 
with a clean image containing photon events with energies between 0.3 and 10 keV. 
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3. Source detection and building of tlie catalogue 



The first step in the detection process consists in the choice of an energy band from which ex- 
tract photons to pass to the detection algorithm. The output of the reduction step is a photon event 
file containing photons with energies in a wide energy range (0.3 - 10 keV). However, it is not 
convenient to use the whole band in order to detect sources. In fact, unless we are dealing with 
nearby massive systems, clusters emit generally very few photons in the hard energy band (with 
energies larger than 4-5 keV), therefore in order to maximize their contrast with respect to the 
X-ray background and enhance the detection probability, it is necessary to work only with the 
softest X-rays. In particular, Scharf (2002) calculates the maximum, or optimal, signal-to-noise 
energy band for galaxy system X-ray emission detected by the Chandra observatory. He finds 
that the "classical" ROSAT-like 0.5-2.0 keV band is close to optimal for clusters in a wide range 
of intracluster medium (ICM) temperature (in particular for T> 2 keV) and redshift z < 1 . Based 
on these considerations, for the detection of sources I choose to consider only photons in the 
energy range 0.5-2.0 keV. 

As a detection algorithm I use the Voronoi Tessellation and Percolation technique (hereafter 
VTP). The VTP method of Ebeling (1993, see also Ebeling & Wiedenmann 1993) is a general 
method for the detection of non Poissonian structure in a distribution of points. By definition, a 
Voronoi tessellation on a two-dimensional distribution of points (called nuclei) is a unique plane 
partition into convex cells, each of them containing one, and only one, nucleus together with the 
set of points which are closer to that nucleus than to any other one (Voronoi 1908). In X-ray 
observations, any nucleus is a detector pixel where at least one photon count occurs. To each of 
these pixels, VTP associates a cell polygon whose sides form the perpendicular bisectors of the 
non-crossing vectors joining the nearest neighbor photon pixels. Since most cells contain only 
one photon, the surface brightness associated with this photon equals the inverse of the product 
of cell area and local exposure time. Multiple counts pixels are rare and associated only to the 
brightest sources. These pixels are flagged and used in the percolation algorithm described be- 
low. 

In Fig. 13 I show the Voronoi tessellation of a typical Chandra ACIS-I photon distribution. The 
sources in this field are apparent to the eye as clusters of small cells. After the implementation 
of the Voronoi tessellation, the cumulative distribution of the inverse areas of the photon cells is 
computed and compared with that expected from a random (Poisson) distribution (Kiang 1966), 
and a cutoff can be determined that defines the global background count for that field. A spatial 
percolation algorithm is then run on the cells with area smaller than a given threshold above the 
background, grouping them according to the excess above the background density and forming 
sources. For a better background estimation the latter two steps are performed iteratively so that, 
as sources are detected, the background estimate is revised and the source groupings are redeter- 
mined. Finally, for each source I obtain a series of parameters: the source position (computed as 
a weighted sum over photons), the detected source count rate, the detected source area, the min- 
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Fig. 4. Voronoi tessellation of a typical Chandra ACIS-I photon distribution. Photons are shown 
as dots and sources are immediately apparent to the eye as clusters of small tessellation cells. 

imal and maxima moments of inertia of the photon distribution, an estimate of the background 
count, and the probability of the source being a statistical fluctuation. Moreover, the fuU set of 
photons associated with each source is stored. 

The VTP procedure is implemented on a pipeline composed by the CIAO task vtpdetect and a C 
code representing an adaptation of the code developed by Ramella et al. (2001) to detect galaxy 
clusters in optical photometric surveys. 

3.1. Source deblending 

As a first step, VTP is run on each field using as a surface brightness threshold the computed 
value of the background surface brightness. However, in general it is not convenient to run VTP 
with a unique surface brightness threshold. In fact, VTP has, by design, a potentially serious 
drawback: it is unable to separate very close pairs of sources, i.e, it inevitably produces "blends" 
as soon as it finds the high-flux zones of two emission regions to be touching one another. In 
order to solve this problem, Scharf et al. (1997) show that a practical solution is to run VTP using 
diff"erent surface brightness thresholds. In this way, it is also possible to reduce the uncertainties 
in source identification resulting from positive background fluctuations that become grouped 
with source photons. I run VTP three times with three diff'erent thresholds (denoted as factors 
of the background surface brightness): 1.0 (the first step), 1.3 and 1.7. Then, I select the sources 
matching the following criteria: 



8 W. Boschin: A deep cluster survey in Chandra archival data. First results 

- Sources at threshold 1.0, which do not deblend at higher thresholds, and with an observed 
count rate greater than 2 x 10"'' counts s"' (corresponding to an observed flux limit of ~ 1 

X 10"'^ erg s"' cm"^); 

- Sources matching the above count rate criteria at threshold 1.3 that were parts of blends at 
1.0 but do not deblend at threshold 1.7; 

- Sources matching the position and count rate criteria at threshold 1 .7 that were parts of blends 
at threshold 1.3. 

This deblending procedure provides a list of candidates that is used in a visual inspection of each 
field. In total, I find 3340 X-ray sources in the 81 selected fields at the threshold 1.0. Raising the 
threshold level, the number of detected sources at thresholds 1.3 and 1.7 is, respectively, 94% 
and 84% of the total number detected at the lowest threshold. 65 sources (2%) are deblended 
with the described procedure. For some individual sources, it was necessary to change by hand 
the detection threshold in order to exclude photons that were clearly positive noise and were in- 
cluded in the objects detected at threshold 1 .0. These spurious photons clearly act to bias the area 
estimate of the source. In general, increasing the threshold removes the noise and typically only 
removes less than 10% - 20% of the source photons. Moreover, these photons will be recovered 
in the flux correction step described below. The threshold used for each source is recorded and 
used in the correction from detected to total flux, as detailed in Sect. 13. 31 

Figs.|5]and|6ldemonstrate the differences between the results for the lowest and highest threshold 
for a portion of an ACIS-I field. Photons identified as belonging to sources by VTP are plotted 
in heavy type. The source numbered 2 in Fig. |5]is a clear blend of two pointlike sources when 
observed at the lowest threshold. In Fig.|6l these two sources are deblended at the higher (1.7) 
threshold. 

It is clear that, in principle, there may be cases of real physical systems containing substructures 
that become deblended. However, a visual inspection of the 81 fields does not reveal any case in 
which deblending had split up what was probably a single extended object into separate compo- 
nents. By visual inspection I also reject a set of 102 sources which appear as border effects or 
artifacts due to strong X-ray sources contained in the selected fields. 

3.2. Separation between pointlike and extended sources 

Once the deblending process is completed, the next step consists in separating between point- 
like and extended sources. Really, the very small value of the PSF radius (less than 0.5 arc- 
sec at the aimpoint), makes this operation quite an easy task when considering on-axis sources. 
Nevertheless, a smart procedure is necessary when considering off-axis sources, expecially where 
the PSF radius is larger than about 10 arcsec. This is expecially true for sources located in chips 
12 and 13 during ACIS-S observations. In order to build such a procedure, assumptions must be 
made about the nature of the sources. In particular, I assume that sources are either intrinsicafly 
extended or pointhke. In modeling a pointhke source I assume the PSF is Gaussian. The VTP al- 
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gorithm returns estimates of the mean surface brightness of the background, (Tbkg, a background 
corrected estimate s of the observed count rate for each detected source, and the area of the source 
above the surface brightness threshold. From the area, I obtain an area equivalent source radius 
defined as ryxp = (Avtp/tt)'^^. Supposing now the source is pointlike, the detected and background 
corrected count rate s can be written as 

s = i(rvTp) = 27!- I o-fSF(r)rdr, (1) 
Jo 

where crpsp(r) = o-o&psfir) is the assumed surface brightness profile for a Gaussian PSF, with 
d-pspif) = exp(-r^/2p^). p is the width of the PSF. Besides Eq. Q, according to Ebeling et al. 
(1996), I need a second equation specifying how close to the background the surface brightness 
of the outermost regions of a VTP detection is: 

0"psf('"vtp) - CTo (5-psf('-vtp) = (/-!) CTbkg. (2) 

Here / is the normalized distance of the lowest surface brightness region to the background 
level. This parameter is returned by VTP for each source. Combining Eq. Q and Eq. it is 
now possible to eliminate cry and obtain an equation in the unknown p which can be solved. I 
stress that at no stage of this procedure the model profile from Eq. Q is actually fitted to the 
observed surface brightness distribution. Although I do assume a specific model (Gaussian in 
this case), it is only its integral properties that enter; thus, even if the ACIS PSF profile is not 
exactly a Gaussian, deviations of the true distribution from the assumed model have much less 
effect on the results than they do for the fitting procedures employed by conventional detection 
algorithms. However, the lack of any radial fitting also entails that the value for the source width 
p determined in this process cannot be expected to be as accurate as what might be obtained in a 
detailed imaging analysis of pointed data. 

Now, the estimated p of the detected sources can be compared with the best-fit PSF p as a function 
of the off-axis angle in order to discriminate between pointlike and extended sources. Of course, 
what happens is that extended sources will have estimated p larger than the expected value for 
a pointlike source. In particular, the lower limits of p for a source to be considered as extended 
are computed by building Monte Carlo simulations of ACIS fields. Sources are simulated with 
pointlike surface brightness, convolved with the instrumental PSF (for a nominal photon energy 
of 1.5 keV) and added to a representative background of Poissonian noise. For this purpose 1 
use the package MARX (e.g. Wise 1997), a set of executables designed to be run in sequence 
to produce a Chandra/ACIS simulation. In particular, I simulate more than 500 pointlike sources 
with different count rates and signal-to-noise ratios SNR (as defined in Sect. |5|i and located at 
different off-axis angles. Figs. 0|8l and |9l show the estimated p for all the simulated pointlike 
sources vs. off-axis angles in the SNR bins >10, 6-10 and <6 respectively. For each figure, the 
solid line is a polynomial model fit to all data points, while the dotted line is a 99.9% confidence 
level line derived from the p residuals of the PSF best-fit model. For each source, the errorbars 
in the p estimate are Icr errors derived from the fit. A source is considered as extended if its p is 
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larger than the PSF p at the 99.9% confidence level in the corresponding SNR range, 
from the 81 selected fields 1 classify 51 extended sources. 



In total, 



3.3. Source count rates correction for missing flux 

The presence of an X-ray background always limits the emission directly detectable by any 
source detection algorithm to some fraction of the total source flux. Therefore, the raw VTP count 
rates of the detected sources have to be corrected for low surface brightness emission in the far 
wings of the source that has not been directly detected. This is particularly important for extended 
sources. In order to perform this correction I have to assume again a surface brightness profile for 
the detected sources. In particular, for the extended sources I adopt the King's approximation to 
the density profile of an isothermal sphere (King 1962) and assume a source profile of the form 
(Cavaliere & Fusco Femiano 1976) 



where (TY_(r) is the projected surface brightness as a function of radius, ctq is the central sur- 
face brightness and rc is the core radius. I also assume /3 - 2/3. This value of y6 is a good 
approximation for nearby clusters (Jones & Forman 1984). At medium-high redshifts, clusters 
are characterized by values of /3 in the range 0.5-0.8 (e.g. Arnaud et al. 2002; Holden et al. 2002; 
Pointecouteau et al. 2002), thus /3 - 2/3 is intermediate between these values. Now, from the 
knowledge of background surface brightness and of the VTP source characteristics I can numer- 
ically compute (To and rc (see Ebeling et al. 1996 for the details of this procedure), and the true 
total source count rate can be determined from 



In order to test this flux correction technique on the extended sources found in ACIS fields 
I perform Monte Carlo simulations of ACIS fields. By using again the package MARX I 
simulate more than 500 extended sources with a King profile surface brightness convolved with 
the instrumental PSF. Then, the simulated sources are added to a representative background 
of Poissonian noise. An ensemble of sources are made for different sets of typical intrinsic 
parameters (extent, flux, off-axis angle) to determine the expected scatter in the VTP estimates. 
In Figs. [TO] ^2 El El results are presented across a large range of source parameters, from 
extents of ~ 2 to 80 arcsec and effective fluxes of ~ 2 x 10"'^ to ~ 1 x lO"'"* erg s"' cm"^. In 
all plots, the ratio of the raw, detected count rates to the true count rates and the ratio of the 
corrected (see above) count rates to the true count rate are plotted on the y-axes. The effective 
VTP-detected S/N of a source is a good indicator of the reliability of any flux estimate. As 
expected, the efficiency of this technique to recover the true count rate of the sources depends 
significatively on their SNR. Since the median signal-to-noise ratio of the real X-ray source 
detection is 8, I use this to divide the simulation results for presentation. The effect of noise 
is apparent in the larger scatter seen in flux estimates for the low S/N (<8) sources plotted in 



,2-1-3/3+1/2 



(3) 




(4) 



W. Boschin: A deep cluster survey in Chandra archival data. First results 1 1 

Figs.ll2landll3lwhen compared to those of the high S/N objects in Figs.[TO|and[n For sources 
of medium to small extent (i.e. < 40 arcsec) with S/N greater than 8, 1 can recover the true flux 
to whitin 10% over the flux range shown here. For sources of larger extent (of which I might 
expect to see very few in this survey), the fluxes are systematically underestimated, as expected, 
since most of the flux now lies below the surface brightness limits. For the low S/N sources (S/N 
< 8), the same general trends are present but are dominated by the scatter caused by noise. For 
small extents (< 40 arcsec) I can still recover the true flux of the faintest objects to within 10% - 
20%. 

3.4. Computation of source fluxes 

The ACIS field experiences non uniform exposure. In ACIS-I observations, for instance, it can 
vary by as much as 20% - 25% from the field center to 10 arcmin radius off-axis. Therefore, 
exposure maps are built for each pointing (with 4x4 arcsec resolution pixels) using standard 
CIAO threads in order to take into account this effect when computing source count rates. 
After the count rate correction procedure, source count rates are converted into fluxes multiplying 
them by a factor which depends on the global characteristics of the telescope-detector system, 
on the value of the line-of-sight neutral hydrogen column density, on the energy spectrum of 
the source and on its redshift. I compute this conversion factor by assuming that the extended 
sources exibit a thermal Raymond-Smith (1977) spectrum with temperature kT - 5 keV and 
solar abundance ratio 0.4. For a given hydrogen column density, this conversion factor is accurate 
to within 8% for an ICM temperature in the range 1-20 keV and a solar abundance ratio in the 
range 0.2-1. The hydrogen column density along the line of sight of the sources is taken from 
Dickey & Lockman (1990). Redshifts are not known for most of the candidate clusters at this 
stage of the work, therefore the fluxes listed in the cluster catalog (see next section) should be 
considered as provisional. 

Finally, I perform a first optical follow-up of the extended sources by analyzing the corresponding 
R band fields in the Second STScI Digitized Sky Survey (DSS-II, e.g. McLean et al. 2000). These 
optical fields are sufficiently deep to allow me to identify low redshift X-ray bright galaxies 
(typically spiral galaxies) that appear as extended sources in the X-ray image. Of the 5 1 detected 
extended sources, 15 are identified as single galaxies. Therefore, these sources are excluded from 
the final catalogue of candidate clusters (as an example see Fig. ll4> . 

4. The catalogue of candidate clusters 

This survey currently covers a total sky area of 5.55 deg^. Using the VTP method, the sample of 
candidate clusters contains 36 objects with detected flux greater than 1 x 10 erg s ' cm"^. In 
table 1 I list the detected candidates in growing order of right ascension. The object number is 
given in Col. 1 . The equatorial coordinates (J2000.0) of the X-ray centroid are fisted in Cols. 2 
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and 3 in degrees with four decimals (typical centroid uncertainties are of 1-2 arcsec). The total 
unabsorbed X-ray flux in c.g.s. units (erg s"' cm"^) and its 1 cr uncertainty are listed in Cols. 4 
and 5. The estimated characteristic radius r„p is given in Col. 6. Col. 7 lists the signal-to-noise 
ratio of the sources (as defined in Sect.|3- Finally, Col. 8 contains notes for individual clusters. 
In this column, I also list the identifications from the literature. Only five candidates seem as- 
sociated in DSS-II images to galaxy concentrations with respect to the background/foreground 
galaxy distribution (see Fig. ll5t and are probably low-redshift systems. Three candidates are re- 
discoveries of systems at low-intermediate redshift detected in previous X-ray surveys. In most 
cases no optical counterpart (galaxy, group, etc.) is seen. This suggests that most candidate clus- 
ters are probably at redshifts z > 0.2 - 0.3, thus invisible in the relatively shallow fields of the 
DSS-II. 

5. Sky coverage of the survey 

As already pointed out in the introduction, the aim of a new X-ray selected cluster catalogue is 
twofold: providing a set of targets for optical/NIR follow-up observations of single interesting 
objects and to address global statistical properties of the cluster population useful for cosmolog- 
ical studies. This second purpose requires the evalution of the statistical properties of the survey 
sample (e.g. the logA^-logS' and the luminosity functions). To do this, the knowledge of the ef- 
fective sky coverage, i.e. the survey area in which all sources above a given flux are detected, is 
necessary. 

Since this survey uses pointed data, each field has a different detection sensitivity to objects of 
a given extent and flux, according to exposure and background. In order to estimate this, VTP's 
detection sensitivity can be measured parametrizing a detection criterion based on a definition 
of the source signal-to-noise ratio. In particular, from simulations and real data I find that the 
criterion for source detection is approximable as 



where nyxp is the total number of photons (source and background) that lie within the equivalent 
source radius r„p - {Ayjp/ny^^, where Ayxp is, as previously defined, the source area within which 
the surface brightness exceeds the VTP surface brightness threshold (Tvtp, which is defined by the 
background (crt,kg) and threshold (/) as o-yjp = f crhkg if > 1-0). Assuming again for the sources 
a King profile as in Eq. 0, it is possible to compute the SNR for sources with given fluxes and 
extents in fields with a given exposure time and background surface brightness (see Scharf et 
al. 1997 for the details of this procedure). Exposure and background information are available 
for each of the 8 1 fields used in this survey, therefore I use them to integrate the sky coverage 
over the area of each ACIS field (including also the PSF variation with off-axis angle) adopting 
the criterium introduced above to determine the detection sensitivity. The final result shown in 
Fig. El is the combined sky coverage of all fields used in the survey to the hmiting threshold 



SNR = 




n. 



> 3, 



(5) 
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(/ = 1.0)'. In this figure I also plot two curves representing the loci, with varying redshift, of 
galaxy groups (Lx = 1 x 10"*^ erg s"', = 100 kpc, dashed line) and galaxy clusters (Lx = 1 
X 10'*'' erg s = 150 kpc, solid line). Redshifts therefore increase right-to-left along these 
curves. 

6. Expected redshift distribution of the detected clusters 

An important issue to discuss when compiling a cluster survey is the estimate of the cluster 
redshift distribution. This is the comoving volume per unit redshift and solid angle (dV/dzdO.) 
times the comoving density of clusters ncom with luminosities above the survey detection limit 
L'™. I write it as 



where (f)(Lx) is the cluster X-ray luminosity function (XLF), i.e. the function describing the num- 
ber of clusters per unit volume with luminosity in the range Lx ^ Lx+ dLx- Now, dV com I dzdO., 
the comoving volume per unit solid angle and redshift, is easily computable once a given cos- 
mology is fixed, l)™ is a function of redshift and depends on the assumed cosmological model, 
on the shape of the cluster surface brightness profile and on the exposure times of the observa- 
tions. I compute L}^ at different redshifts assuming, as in the previous section, a King profile 
for clusters, a typical exposure time of 40 ks (which is close to the median value of the exposure 
times of the selected pointings) and SNR > 3 as a criterium for cluster detection. For simplic- 
ity I also assume a cluster core radius of 150 kpc, a quite common value (e.g. Ota & Mitsuda 
2002). Finally I have to integrate the cluster XLF. Following Rosati et al. (2000), I assume it is 
an evolving Schechter function: 



with = ^xo(^ ^) ' where A and B are two constants parametrizing the possible evolution 
of the XLF with redshift. For the local XLF, I assume the one derived from the BCS sample 
(Ebeling et al. 1997). The results are shown in Fig.^] The three curves represent the expected 
cluster redshift distribution computed for three different couples of the evolutionary parameters: 
A = -0.5 and B = -1 (dotted hne), A = -1 and B = -2 (dashed hne), and A = B = 
(no evolution, solid line). These curves show that the expected cluster redshift distribution has a 
maximum at intermediate redshift, at z ~ 0.3 - 0.4. However, a significant fraction of the clusters 
(about 30-40% for the two evolutionary models, higher for the no evolution case) are expected 
to have redshifts higher than 0.6-0.7, with the possibility of having clusters at z > 1. The de- 
tails of this result are dependent on the assumptions done above on cluster profiles. However, I 
find that changing the cluster profile parameters does not affect the quahtative behaviour of the 
cluster redshift distribution. Of course, at high redshifts I expect to discover only relatively rich 
' For calculations in this and in the following sections cosmic K-corrections from Jones et al. 1998 are 
used. 



'X 




(6) 



<^(Lx) = <^o (1 + zt L^" exp(-Lx/L*^), 



(7) 
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systems, while at low and intermediate redshift the contribution of small galaxy groups becomes 
significant. This is evident from Fig.[^ where the fractional sky coverage is shown as a function 
of redshift for three classes of objects. I have complete sky coverage for clusters with Lx = 1 
X 10"*^ erg s"' and core radius - 150 kpc out to a redshift z ~ 0.8, and 50% coverage at a 
redshift z ~ 1.3-1.4. Instead, 1 have a good coverage for groups (or faint, small clusters) of 
Lx = 1 X 10"*^ erg s"' out to z ~ 0.5 (50% coverage). Galaxies have a > 50% coverage out to 
z ~ 0.2 - 0.3. Therefore, galaxies observed in Chandra fields are generally low redshift objects. 
They can be quite easily recognized in the fields of the Digitized Sky Survey and, eventually, 
excluded from the catalogue of extended sources, as explained in Sect. 13. 41 
The knowledge of the sky coverage also allows to compute the volume that the survey probes 
above a given redshift z, for a given X-ray luminosity, e.g. the characteristic luminosity LJ^ cor- 
responding to the "knee" of the XLF. Rosati et al. (2002) compute the search volumes of four 
cluster surveys: EMSS (Gioia et al. 1990), 160 deg^ survey (VikhHnin et al. 1998), NEP survey 
(Henry et al. 2001) and RDCS (Rosati et al. 1998). In Fig.[T9ll plot the volume probed by these 
four surveys as a function of redshift together with the volume probed by this survey. Since these 
surveys cover different solid angles at varying fluxes, they probe different volumes at increasing 
redshift and therefore different ranges in X-ray luminosities at varying redshifts. From the plot 
we see that the EMSS is very sensitive to the most luminous systems, but is able to find only 
a few clusters at high redshifts due to its bright flux limit. NEP, 160 deg^ and RDCS are much 
deeper surveys and discover many clusters at high redshifts. In particular, the RDCS pushes the 
search to the faintest fluxes, providing sensitivity to the highest redshift systems with character- 
istic luminosity even beyond z = 1. Instead, the NEP and the 160 deg^ surveys do not reach 
the flux limit of RDCS and probe the bright end of the XLF out to z ^ 1 . This survey (the thick 
sohd line in the figure) probes a small total volume with respect to the other surveys bacause of 
the still limited surveyed solid angle. Nevertheless, it is superior to NEP and 160 deg^ surveys in 
probing the Universe at z > 0.8 and competitive with RDCS for redshifts z > 1 . 

7. log A^-log S relation for the candidate clusters 

The computed sky coverage of the survey is useful to derive the statistical properties of the 
cluster population. In particular, it is interesting to calculate the cumulative number density of 
the candidate clusters as a function of flux (the so-called log A^-log S relation) and compare it 
with previous studies. In order to perform this calculation, each cluster is added to the cumulative 
distribution with the weight equal to the inverse solid angle corresponding to its measured flux 
and angular core radius. The derivative cumulative logA^-logS' function is shown in Fig. 1201 
Counts are shown as filled circles. The error bars represent the Icr uncertainties in the number 
of detected clusters. Shown on abscissa is the total source flux, where "total flux" refers to the 
flux extrapolated to surface brightnesses below the detection limit (see Sect, \3.3l . I also show 
the cluster counts derived from other two deep surveys: The RDCS (Rosati et al. 1995, crosses) 
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and an ultra deep UK ROSAT survey (McHardy et al. 1998, asterisk). The logA^-log5 relation 
derived from this survey spans a flux range from 7 x 10"'^ erg s"' cm"^ to 3 x 10"'"^ erg s"' 
cm"^. The lower flux limit is set at a level in which the sky coverage is about 1 deg^ and all 
clusters are expected to be detected as extended sources since they have core radii larger than 
about 15 arcsec, based on a range of typical sizes and luminosities (see Fig. \161 . Therefore, I 
expect this sample be complete to this flux level. At the bright end, the results are in agreement 
within errors with the deep sample from Rosati et al. (1995). At the faintest end, the extrapolation 
of the logA^-logS' agrees with the results of McHardy et al. (1998), who identified most of the 
X-ray sources, regardless of extent, in their ultradeep ROSAT survey. 

8. Summary 

I present a catalogue of 36 candidate clusters detected as extended sources in 8 1 Chandra ACIS 
observations covering an angle of 5.55 deg^. To detect these sources, I use a detection algorithm 
based on the Voronoi Tessellation and Percolation technique. This is an efficient method in order 
to detect low surface brightness extended sources in a shape independent way. Despite the small 
solid angle covered so far, the survey probes a volume at high redshift (z > 0.8) comparable 
or larger than previous Einstein and ROSAT-based surveys. Moreover, I expect a significative 
fraction of the candidates are medium-high redshift (z > 0.6 - 0.7) systems. This makes these 
candidate clusters interesting targets for optical/NIR follow-up observations. In particular, the 
confirmed distant clusters will be useful, for instance, to study the properties of the cluster galaxy 
population at high redshift (e.g. Gioia et al. 1999; Stanford et al. 2002). I also present a prelim- 
inary logA^-log5 relation derived from the catalogue. This relation shows a general agreement 
with some deep ROSAT-PSPC based cluster surveys. As soon as new data will be available in the 
Chandra archive, they will be reduced in order to increase the sky area and the volume probed by 
the survey and enrich the catalogue compiled so far. 
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Fig. 5. Source photons and sources identified by VTP at threshold 1.0 in a portion of an ACIS- 
I field. Thick points indicate VTP source photons, sources are labeled numerically. Note that 
source 2 is a clear blend. 
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Fig. 6. Source photons and sources identified by VTP at threshold 1 .7 in the field shown in Fig.|5l 
Note that source 2 at threshold 1.0 has now been resolved into sources 2 and 3. 
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Fig. 7. Estimated p and 1 cr errors for a set of pointlike simulated sources with SNR>10 vs. 
ofF-axis angle. The solid line is a polynomial fit to all data points. The dotted line is a 99.9% 
confidence level line derived from the p residuals of the PSF best-fit model. 
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Fig. 8. Estimated p and 1 cr errors for a set of pointlike simulated sources with 6<SNR<10 vs. 
oflF-axis angle. The solid line is a polynomial fit to all data points. The dotted Une is a 99.9% 
confidence level Une derived from the p residuals of the PSF best-fit model. 
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Fig. 9. Estimated p and 1 cr errors for a set of pointlike simulated sources with SNR<6 vs. off-axis 
angle. The solid line is a polynomial fit to all data points. The dotted line is a 99.9% confidence 
level line derived from the p residuals of the PSF best-fit model. 
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Fig. 10. Mean ratios of detected (filled circles) and corrected (crosses) count rates to the true 
count rate of high SNR (>8) simulated sources plotted against their extent (arcsec). 
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Fig. 11. Mean ratios of detected (filled circles) and corrected (crosses) count rates to the true 
count rate of high SNR (>8) simulated sources plotted against their intrinsic flux (c.g.s. units). 
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Fig. 12. Mean ratios of detected (filled circles) and corrected (crosses) count rates to the true 
count rate of low SNR (<8) simulated sources plotted against their extent (arcsec). 
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Fig. 13. Mean ratios of detected (filled circles) and corrected (crosses) count rates to the true 
count rate of low SNR (<8) simulated sources plotted against their intrinsic flux (c.g.s. units). 
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Fig. 14. This figure represents an R-band DSS-II image of an ACIS-I field. The ellipses repre- 
sent the sources detected by VTP in the corresponding X-ray image, two of which are clearly 
extended. The image shows they correspond to relatively nearby galaxies. Therefore, these ex- 
tended sources are not listed in the catalogue of candidate clusters. 
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Fig. 15. X-ray and DSS-II R band images of the ACIS-S field containing the candidate cluster 
12. Ellipses indicate sources detected by VTP in the X-ray image. The largest one indicates the 
cluster. Note, in the optical image, the presence of a concentration of faint galaxies in corre- 
spondence with the position of the candidate cluster. It suggests the existence of a low redshift 
system. 
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Table 1. The cluster catalogue 



Number 


R.A. (J2000) 


Decl. (J2000) 


Fx 




<5Fx 




rvTp 


SNR 


Note 




degrees 


degrees 


c.g.s. units 


c.g.s. units 


arcsec 






(1) 


(2) 


(3) 


(4) 




(5) 




(6) 


(7) 


(8) 


1 


7.7555 


+26.4441 


2.1x10- 


14 


0.3 X 10- 


14 


26.2 


7.73 




2 


24.1793 


+20.8143 


1.0 X 10- 


14 


0.2 X 10- 


14 


29.6 


6.85 


P 


3 


44.3881 


-23.4017 


0.3x 10- 


14 


O.lx 10- 


14 


17.7 


2.91 




4 


46.4686 


+03.8439 


0.8x10- 


14 


0.3x10- 


14 


21.3 


3.82 




5 


54.4602 


-23.0730 


2.4x10- 


14 


0.3 X 10- 


14 


27.2 


9.52 




6 


73.5050 


-03.1444 


2.0x10- 


14 


0.3 X 10- 


14 


27.3 


6.42 




7 


73.6638 


-03.1277 


0.7x10- 


14 


0.2x10- 


14 


16.6 


3.84 




8 


73.7697 


- 10.2003 


3.2x10- 


14 


0.3 X 10- 


14 


49.3 


10.95 




9 


73.7895 


-10.1752 


1.4x10- 


14 


0.2x10- 


14 


37.5 


6.71 




10 


85.6665 


-40.9155 


0.5 X 10- 


14 


0.1x10- 


14 


14.7 


6.63 


P 


11 


125.0098 


+63.7611 


1.7x 10- 


14 


0.3 X 10- 


14 


30.2 


5.79 




12 


133.0950 


+51.3794 


2.1 X 10- 


14 


0.2 x 10- 


14 


36.4 


13.45 




13 


133.2558 


+51.2218 


1.2 X 10- 


14 


O.lx 10- 


14 


26.1 


11.33 


B 


14 


137.5350 


+54.3149 


1.4 X 10- 


14 


O.lx 10- 


14 


27.4 


15.45 




15 


142.7633 


+79.2219 


3.0x10- 


13 


0.1x10- 


13 


75.5 


29.45 




16 


148.6100 


+68.9119 


0.7x10- 


14 


0.1x10- 


14 


26.5 


6.37 


B 


17 


158.8064 


+57.8383 


1.0x10- 


14 


0.1x10- 


14 


29.3 


7.74 




18 


158.8564 


+57.8472 


1.2x10- 


14 


0.1x10- 


14 


29.1 


8.69 




19 


164.0525 


-03.5845 


5.0x10- 


14 


0.2x10- 


14 


33.7 


19.27 


P 


20 


169.3591 


+07.7274 


2.6x10- 


14 


0.3 X 10- 


14 


36.9 


8.84 


1 


21 


169.3746 


+07.7716 


0.9x10- 


14 


0.2x10- 


14 


26.7 


4.77 


2 


22 


175.0341 


-26.5263 


1.5 X 10- 


14 


0.2x10- 


14 


32.0 


8.23 




23 


199.2258 


+29.2394 


2.4 X 10- 


15 


0.5 X 10- 


15 


19.5 


6.05 




24 


199.3051 


+29.3089 


1.2 X 10- 


15 


0.4 X 10- 


15 


14.9 


4.14 




25 


204.4194 


+29.4313 


2.0 X 10- 


15 


0.8 X 10- 


15 


17.3 


3.14 




26 


209.8918 


+62.3164 


3.1x 10- 


13 


0.6 X 10- 


14 


69.4 


49.72 


3 


27 


218.1304 


-01.1376 


1.4x 10- 


14 


0.3x10- 


14 


26.4 


5.82 




28 


233.5057 


+23.6009 


1.1x10- 


14 


0.1x10- 


14 


19.7 


8.21 




29 


234.9979 


-03.0390 


5.7x10- 


14 


0.3 X 10- 


14 


51.1 


19.89 


B 


30 


237.3827 


+21.5499 


3.4x10- 


14 


0.3 X 10- 


14 


33.4 


14.15 




31 


245.9976 


+26.6492 


0.5 X 10- 


14 


0.1x10- 


14 


14.5 


4.50 




32 


249.1558 


+41.1306 


0.8 X 10- 


14 


0.1x10- 


14 


25.6 


9.85 




33 


255.7827 


+51.5614 


0.9x10- 


14 


0.2x10- 


14 


29.3 


5.08 


B 


34 


310.9874 


+77.1801 


0.7x 10- 


14 


0.2 X 10- 


14 


24.1 


3.46 




35 


333.2226 


-22.1602 


l.lx 10- 


14 


0.2 X 10- 


14 


23.7 


6.03 




36 


345.7229 


+08.8118 


1.3 X 10- 


15 


0.4x 10" 


15 


12.4 


3.81 





NOTES - (1) z=0.4 cluster (N. 97 Vikhlinin et al. 1998), (2) z=0.16 cluster (RIXOS F258_ 101, Mason et al. 
2000), (3) cluster MS 1358.4+6245 at z = 0.328, (P) flux uncertain due to the presence of a nearby pointhke 
source, (B) source very close to the border of the X-ray image. 
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Fig. 16. Sky coverage of the 81 fields used in this survey. Sky coverage is plotted as a function of 
intrinsic flux and intrinsic, projected core radius (assuming King profile). Contours are drawn at a 
percentage of the total survey area. From the left, the first contour is at a level of 1 %, subsequent 
solid contours are 10%, 20%,..., 90%, 95% and 100%. The two curves at the right-hand side of 
the figure are the loci, with redshift, of galaxy groups (Lx - 1 x 10'*^ erg s"', rc - 100 kpc, 
dashed Une) and galaxy clusters (Lx = 1 x 10^"* erg s - 150 kpc, solid line). Redshifts for 
these systems increase right to left. The points corresponding to some specific redshift values are 
marked and labeled on the two curves. The vertical dashed line corresponds to the flux value 7 x 
10"'^ erg s"' cm"^ (see Sect.Q. 
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Fig. 17. Expected redshift distribution of the clusters detected in this survey in the redshift range 
< z < 1.5. Three curves related to three different couples of the XLF evolutionary parameters 
(see text) are drawn. 
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Fig. 18. Sky coverage as a function of redshift for three classes of objects: elliptical galaxies 
(with Lx = 1 X 10"*^ erg s"' and effective core radius = 30 kpc) in dashed line, groups (Lx = 
1 X 10"*^ erg s \ Tc - 100 kpc) in dotted line, and clusters (Lx = 1 X 10"^"* erg s"\ rc = 150 kpc) 
in soUd line. 
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Fig. 19. Search volumes of five X-ray cluster surveys, V(> z), for a cluster of given X-ray lumi- 
nosity (Lx - 3 X lO^'^/Zjg erg s ' ^ L^). The thick solid line corresponds to the volume probed 
by this survey (adapted from Rosati et al. 2002). 
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Fig. 20. Cumulative logA^-log5 function of the candidate clusters. Counts are shown as filled 
circles. The error bars represent the Icr uncertainties in the number of detected clusters. Data 
from Rosati et al. (1995, crosses) and from McHardy et al. (1998, asterisk) are also shown for 
comparison. 



